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Abstract 

We explore possibility of tuning photonic crystal properties via order-disorder transition. We 
fabricated a photonic bandgap material consisting of a three-dimensional array of conducting mag- 
netizable spheres. The spheres self-assemble into ordered state under external magnetic field, in 
such a way that the crystalline order can be continuously controlled. We study mm-wave trans- 
mission through the array as a function of magnetic field, i.e. for different degrees of order. This 
was done for the regular crystal, as well for the crystal with the planar defect which demonstrates 
resonance transmission at a certain frequency. We observe that in the ordered, "crystalline" state 
there is a well-defined stopband, while in the completely disordered, glassy or "amorphous" state, 
the stopband nearly disappears. We relate the disappearance of the stopband in the disordered 
state to the fluctuations in the particle area density. We develop a model which predicts how these 
fluctuations depend on magnetic field and how they affect electrodynamic properties of the whole 
sample. The model describes our results fairly well. 
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INTRODUCTION 



Photonic bandgap materials are ordered arrays of scatterers that do not allow electro- 
magnetic wave propagation in some frequency ranges named gaps pQ. The important chal- 
lenge in the field of photonic crystals is tunability, in other words, possibility to control the 
depth, width and position of these frequency gaps. This can be achieved by various means 
such as liquid crystal infiltration [2], temperature elastic strain [I], and magnetic field 
El [7J IH [9j [TOj [HI [121 1131 E] • We explore here a novel route to achieve tunable photonic 
bandgap materials: magnetic-field-induced order- disorder transition. 

Since the gaps in photonic bandgap materials arise from their periodically-ordered struc- 
ture, we would apriori expect that disorder destroys the gaps. However, it is well-known 
that disorder by itself may result in wave localization, i.e. it may create the gap [15J. Hence 
the effect of disorder on photonic bandgap materials may be multidirectional and it comes 
at no surprise that it has been studied so intensively during last years. Numerous analytical 
and numerical studies showed that the gaps in photonic crystals are robust [TBI El HE], i-e. 
they do not disappear under weak and moderate disorder. More specifically, it was found 
that under weak disorder the gap edges become smeared [191 EDI ETJ E2J E31 EH [25], the 
higher-order gaps become more shallow [26], while the depth of the first stopband remains 
almost unaffected [T71 HH [20J (2TJ [221 1211 [26]. Under moderate disorder the first stopband 
becomes more shallow [121 122] • Under strong disorder there appear localized states in the 
gaps, which are characterized by enhanced transmission [2"01 f2T\ |2"1]. Theoretical predictions 
indicate that the effect of disorder is most pronounced at the gap edges and is minimal in 
the passbands, at the frequencies corresponding to the so-called Bragg remnant or antigap 
[261 E3 EH [29] . Various kinds of disorder (positional disorder, size disorder, fluctuations in 
refraction index, etc.) produce almost the same qualitative effect. The only exception is 
sliding disorder (or stacking faults) which can increase the depth of the gap [221 EH]- 

There are only a few experimental studies of the wave propagation in photonic bandgap 
materials with disorder. Ref. [21J studied microwave propagation through the two- 
dimensional arrays of rods and found bandgap smearing upon increasing positional dis- 
order, while Ref. [30] studied microwave transmission through random assemblies of spheres 
and observed wave localization and stopband resulting from disorder. Introduction of small 
controlled disorder (in fact, " quasicrystalline" order) into two-dimensional photonic bandgap 
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material allows to achieve a full photonic stopband in the near- infrared [3T| [32] . 

In this work we explore possibility of tuning of photonic bandgap materials by continuous 
variation of the degree of disorder. This is done in the mm-wave range with a model system 
of metallic magnetizable spheres whose lateral position is controlled by external magnetic 
field. This resembles the polymer-dispersed liquid crystals where fluctuations in the optical 
birefringerence are controlled by electric field [33]. A similar idea (although based on particle 
reorientation rather than motion in magnetic field) has been suggested with respect to 
photonic bandgap materials [8]. 



EXPERIMENTAL SETUP 



Photonic crystal with tunable disorder 

We build our photonic crystal from the 2 mm diameter steel spheres. In the presence 
of magnetic field each sphere acquires magnetic moment proportional to the field, M = 
3xVH/2 where H is external magnetic field, V is the volume of the sphere, and \ is the 
magnetic susceptibility. We used SQUID magnetometer to measure particle magnetization 
in the fields up to 500 Oe and found coercive force of 5 Oe and volume susceptibility x=0.15. 

We put 397 spheres into a 0.67 mm thick plexiglas container (Figs.lj2) and mounted 



several such containers in the stack. The interlayer spacing is 3.5-4.5 mm, while the nearest- 
neighbor distance in a single layer is 3 mm. The stack is mounted inside the Helmholtz coils. 
In the presence of external magnetic field the spheres become magnetized in the direction 
perpendicular to the layers, in such a way that the spheres in each layer repel each other. 
The out-of-plane attraction between the spheres is small compared to the in-plane repulsion. 

The static particle configuration in each layer is determined by the interplay between (i) 
magnetic forces and (ii) static friction forces between the spheres and the substrate. Under 
strong magnetic field, magnetic forces exceed friction forces and the particles self-assemble 
into hexagonally-ordered crystalline lattice whose orientation is determined by the container 
shape. Since we use the hexagonally-shaped container, and the total number of particles 
corresponds to the perfect hexagonal packing, the resulting lattice is almost perfect at strong 
magnetic field - Figjl] At the intermediate magnetic field, the friction forces deform this 
crystalline lattice, in such a way that it splits on several crystalline grains separated by grain 
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boundaries. In the weak magnetic field, the friction forces dominate, the particles are in the 
disordered (amorphous or glassy) state and may touch one another. 

To characterize disorder we took images of the particle configuration using a CCD camera 
and then performed digital Fourier-transform of the images. In the strong field we observe 
sharp Bragg peaks (FigjTJ right panel) indicating on crystalline order; at the intermediate 
field the intensity of these peaks decreases and the diffuse background appears (Fig. lb); 
in the weak magnetic field the Bragg peaks disappear, indicating on the absence of the 
crystalline order (Fig. [T|. 

The images of our arrays (FigJI]) are very similar to the snapshot of the particle config- 
uration in the system of interacting magnetic dipoles at finite temperature [35]. However, 
instead of temperature, the disordering agent in our case is static friction between the par- 
ticles and the substrate. 

Mm-wave transmission measurements 

We measured mm-wave transmission through our device at normal incidence and in the 
range of 20-50 GHz. We used a HP850C Vector Network Analyzer as a source and two 
standard gain horn antennae to which we attached home-made collimating teflon lenses 
(Figj2|. The inner curvature of the lens is 15 cm, the outer curvature is 30 cm, the distance 
between antennae is 28 cm, the beam diameter is 6 cm. The height of the device is 4-5 
cm. To prevent diffraction at the edges of the stack, we put two apertures above and below 
the stack. The aperture size is slightly smaller than the container diameter. The receiving 
antenna accepts radiation in the angle of 25°, hence it measures forward transmission and 
small angle scattering as well. The noise floor of our setup is -50 dB (measured by replacing 
a sample with a conducting foil of the same diameter). Calibration was performed without 
sample. The size of the particles, the distance between them and the distance between 
the layers was carefully chosen in such a way as to achieve a wide and deep stopband in 
the accessible frequency range. This was done by computer simulations using ANSOFT 
software. 

We studied transmission through the stack of 6-10 layers as a function of magnetic field, 
i.e. at different degrees of in-plane disorder. The measurements were performed as follows. 
We applied a certain magnetic field to initially disordered array and either gently vibrated 
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it or sent a short current pulse ("magnetic stirring") to achieve equilibration. Then we 
measured mm-wave transmission through the stack. At the next step we switched off the 
magnetic field and vibrated the stack once again to recover initial, disordered state. Then 
we repeated the measurement for another value of magnetic field and so on. While there 
exist many particular configurations corresponding to a certain value of the external mag- 
netic field, the mm-wave transmission through these metastable configurations differ only in 
minor details. Therefore, we performed 3-10 vector transmission measurements at the same 
value of magnetic field and for different particle configurations; and vectorially averaged the 
results. We also checked experimentally that the transmission through the stack is almost 
independent of the polarization of the incident wave. A small lateral shift of one of the 
layers with respect to another (sliding disorder) was also found to be unimportant. 

EXPERIMENTAL RESULTS 

Magnetic field dependence of the mm-wave transmission 

Figure [3] shows mm-wave transmission through the six-layer stack at two extreme values 
of magnetic field: one corresponding to completely ordered state, and another corresponding 
to completely disordered state. In the ordered state there is a well-defined stopband at 24-44 
GHz. For the stack with only six layers, the transmission in the stopband is very small, T = 
35 dB. Reflectivity in the stopband is close to unity (not shown here) [31]. Transmission in 
the passband can be as high as -1-2 dB, indicating on negligible absorption losses. This is 
not surprising since the losses are determined by the ratio of the skin-depth 5 in the material 
of the sphere to its radius r. For the mm-wave frequencies this ratio is so small ~ 10~ 3 , that 
the spheres can be considered ideally conducting, i.e. lossless. The spikes in the stopband 
in the Figj3] do not arise from the noise of the measurement system. They are related to 
the finite number of layers and to the resonant transmission. The spikes in the magnitude 
of transmission are accompanied with the similar features in the phase (not shown here). 

Upon decreasing magnetic field, the stopband becomes smeared and eventually disap- 
pears. Transmission monotonously decreases with decreasing field at all frequencies. The 
strongest effect of magnetic field is at the stopband edges where it achieves 30 dB; while in- 
side the stopband the magnetic field effect is negligible. Magnetic field effect in the passband 
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is also quite pronounced. To study it in more details we prepared a similar stack but with 
increased interlayer separation. Here, the stopband is shifted to lower frequencies in such a 
way that the frequency range above the first stopband becomes "visible". Figure [4] shows 
the effect of magnetic field on the mm-wave transmission in this crystal. The magnetic field 
effect in the stopband may achieve 15 dB. Note the broad transmission peak between the 
first and the second stopband in the disordered state, in the frequency range corresponding 
to the passband in the ordered state. This peak of enhanced transmission is known as Bragg 
remnant or antigap 12H I2H] ■ 

Figure [5] shows magnetic field dependence of the mm-wave transmission for the sample of 
Fig.4 and at few different frequencies, corresponding to the midgap, gap edge, and antigap. 
Although the overall variation of transmission strongly depends on frequency, the field de- 
pendence is much more the same: strong linear increase below 60 Oe followed by saturation 
at higher field. The only exception is the midgap frequency where the overall change in 
transmission is so small that its functional dependence cannot be reliably determined. 

This can produce a false impression that the mm-wave transmission is totally insensitive 
to magnetic field for the frequencies corresponding to the midgap. However, this is not 
exactly so. Figure [6] shows mm-wave transmission through the stack in which we introduced 
a planar defect, i.e. we displaced two halves of the stack (Figjl]) in such a way that the 
distance between two central layers is increased as compared to other interlayer distances. 
This results in a sharp transmission peak inside the stopband at 34.7 GHz. Upon decreasing 
magnetic field this peak decreases and completely disappears at H = 0. The effect of 
magnetic field on transmission at the resonance frequency amounts to 30 dB. Therefore, 
magnetic field can strongly affect transmission in the stopband provided this transmission 
is resonant. 

Dispersion relations in the ordered and disordered states 

Magnetic field affects not only the magnitude of the transmission (Figs, 
phase as well. This means that the phase and group velocity can be magnetically tuned. 
These velocities are determined from the dispersion relation k(f), where k is the wavevector 
and / is the frequency. To determine k(f) we performed vector transmission measurements 
at fixed field and for the stack with varying number of layers, N. The real and imaginary 



3p|5|6) but the 
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parts of the wavevector were found from the experimental data using the following relations: 

where T is the power transmission coefficient, is the phase shift upon transmission, d z 
is the unit cell period in the direction of propagation. Figure [7] shows our results. In the 
ordered state, the Im(A;) is very small beyond the stopband, as expected in the lossless 
material. This is in contrast to the disordered state where the lm(k) is more or less the 
same at all frequencies. Note the gap in Re(fc) in the ordered state and the absence of the 
gap in the disordered state. The group velocity is v g = dk/df. In the disordered state it 
is fairly frequency-independent, while in the ordered state it demonstrates strong frequency 
dependence, in particular, v g — > at the stopband edges. Therefore, group velocity can be 
continuously tuned by magnetic field. 



MODELING 



Although particle configuration in each layer strongly depends on magnetic field, to un- 
derstand how the changes in configuration affect the mm-wave transmission through the 
whole sample is not an easy task. Indeed, to explain the appearance of the stopband for 
the certain direction of propagation in a three-dimensional photonic bandgap material, it 
is usually enough to represent it as a multilayer with periodicity only in the direction of 
propagation. In this representation the layers are assumed spatially uniform and the details 
of the particle arrangement in each layer are of little importance. Our idea is that any devi- 
ation from the in-plane crystalline order leads to fluctuations in the area density of particles. 
Therefore, the layer can be no more considered as uniform. 

The purpose of our modeling is to find how the in-layer disorder affects the electro- 
magnetic wave propagation perpendicular to the layers. We attribute the effect of lateral 
disorder to the area density fluctuations. The modeling is performed in three steps. First we 
estimate fluctuations of the particle area density as a function of magnetic field. Secondly, 
we calculate refraction index and reflectivity of a single layer of spheres with fluctuating 
area density. Thirdly, we consider wave propagation through the stack of such layers. 
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Magnetic-field-dependent disorder in a layer of magnetizable spheres 



Particle configuration in our array is determined by the interplay between magnetic and 
friction forces. Consider first magnetic forces. The potential energy of induced magnetic 
dipole-dipole repulsion is: 

M,.\l, , . 

where 7*j is the particle position, Mj is the particle magnetic moment. Magnetic forces 
Fmagn = Ylj drive particles into equilibrium hexagonally-ordered state where magnetic 
repulsion energy is minimized and F % magn = 0. The lattice constant of this array is deter- 
mined by the number of particles and the size and shape of container [T2]. Since magnetic 
susceptibility of our spheres is not high, we assume M oc H where H is external field (rather 
than the sum of the external field and the induced field of other particles). 

To achieve equilibrium, the particles should move laterally. Here, static friction forces 
become important. Indeed, to push a hard sphere on flat substrate out of equilibrium requires 
application of some mimimal force, exceeding static friction force, Ff r = \ij r P where P is 
normal force on a particle and /// r is the rolling friction coefficient which depends on the 
radius of the sphere and on the nature of contacting surfaces. In the presence of static friction 
the particles start to move only when F magn > Ff r and come to rest when F magn < Ff r . 
To characterize the interplay between magnetic and friction forces we introduce "magnetic 
length" 21 h that represents a minimal distance between two isolated magnetized spheres on 
substrate with friction. We find it from the relation F magn (2l H ) = Fj r . Equation (|2j yields 

In this work we are interested in the magnetizable particles (M oc H), hence Ih oc H 1 ^ 2 . 
The magnetic length should be compared to the average distance between the particles. The 
latter is conveniently characterized by the Wigner-Seitz radius a = (7rp)~ 1//2 , where p is the 
area density of the particles. Note that Ih is controlled by magnetic field (Eqj3]) while a 
is field- in dependent. In what follows we consider different regimes, defined by the relation 
between Ih and a. 
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Strong magnetic field (crystalline state) 



Here Ijj » a. Magnetic repulsion tends to keep the particles at equal distances. For 
confined array and at high enough field this results in the hexagonally-ordered "crystalline" 
state. The whole array is a single crystallite (with probably few point defects) which is 
deformed by friction forces. In the absence of point defects such as dislocations and discli- 
nations [36, 37J, we estimate the area density fluctuations arising from elastic deformations 
as follows. The condition of static equilibrium of the elastically deformed planar crystalline 
lattice is 

Here is the stress tensor and fi is the density of bulk external forces which in our case are 
friction forces. The magnitude of these forces is Ff r . We assume that their correlation radius 
is the Wigner-Seitz radius, a. By integrating Eqjl] over unit cell, we find \o~ik\ ~ pFf r a. 
The compressional deformations are \uu\ ~ pFf r a/K, while the shear deformations are 
\ u ik\ ~ pFf r a/fi,(i 7^ k). Here, K is the bulk compression modulus and p is the shear 
modulus. For the planar hexagonal lattice of parallel magnetic dipoles 



2.147M 2 T ^ 1n 
/2 = ;K = Wp (5) 



a 5 



Equations [3] and [5] yield fluctuations of the deformations, \uu\ ~ 2.75 x 10 3 (a//#) 4 ; 
\uik\ ~ 10|ixm|. The corresponding density fluctuations are very small: 

^^(f]\ H ^ (6) 



P 2 " \Ih 
Intermediate field (polycrystalline state) 

Here In > a. In this regime, the friction forces between the particles and the substrate 
are strong enough to split the crystalline lattice onto separate grains with well-defined grain 
boundaries. To estimate the grain size R, we follow Larkin-Ovchinnikov treatment of the 
pinned vortex lattice in superconducting films [38] ■ The maximum energy, associated with 
one particle in the lattice pinned by friction force, is ~ Ff r a. Since the friction forces 
are randomly oriented, the average pinning energy for the grain of area ttR 2 , containing 
Z g = (R/a) 2 particles, is ~ FfraZ^ 2 = Ff r R. These pinning forces lead to displacement of 
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the grain boundary by ~ a. The elastic energy associated with this displacement is ~ pa 2 
The excess energy per unit area of such lattice is: 

5E pa 2 - F fr R 



S ttR 2 



(7) 



Since K » p we neglect here the elastic energy due to compressional deformations. Mini- 
mization of Eq|7] yields the equilibrium grain size, R = 2pa 2 /Fj r . Using Eqs 3j5 we find the 
grain size 

74 

R « 22.44 oc H 2 (8) 
cr 

To estimate density fluctuations we note that the lattice in the grains is almost perfect, 
while the particle configuration at grain boundaries is distorted and strong deformations 
Uik occur there (FigjTJ). In other words, << 1 in the grains and Uik ~ 1 at grain 
boundaries. The relative concentration of grain boundary particles is Zg ~ 2a/ R. The 
average fluctuation of the deformations in the area S, containing Zs = pS particles, is 

-j- Z B 0.09/ a \* 2 
Weak magnetic field (amorphous state) 

Here Ih « a. The particle arrangement reminds a glass or amorphous solid consisting 
of impenetrable spheres with the radius Ih- The area density fluctuation in the area S 
containing Zs particles is [39J: 



5p 2 _ AZ 2 _ 1 + / uOS 



(10) 



P 2 Z 2 Z s 

where v is the pair correlation function. Following Ref . jlO] , we assume that the correlation 
function for the dense disordered planar array of impenetrable spheres of radius r is such 
that J udS ~ —Txpr 2 . We substitute r by Ih, p by 1/na 2 , introduce these values into EgJlOj 
and find 



p 2 z. 



s 



(11) 



Since Ih oc H 1 ^ 2 , the Eq|ll] yields density fluctuations linearly decreasing with H. Note, 
that if Ih < r, where r is the radius of the sphere, the magnetic forces are too small as 
compared to friction forces, and magnetic field is inoperative. 
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Electrodynamic parameters of a planar array of conducting spheres 



We consider first an ordered planar array of ideally conducting spheres of radius r, such 
that r < X/2tt, where A is the wavelength of the incident wave (such spheres can be con- 
sidered as Rayleigh scatterers). The average nearest-neighbor distance is smaller than A/2. 
When the planar electromagnetic wave is incident on such a layer, it is not absorbed but 
scattered. The scattering occurs mostly in the backward and forward directions, while the 
scattering in oblique directions is strongly suppressed. Therefore, for the normally incident 
plane wave this layer can be considered as a uniform non-absorbing medium with thickness 
d = 2r, effective refraction index n and admittance Y [UJ. For dilute arrays of scatterers 
the Refs. pEJ S3] yield 

n^- a ^m (12) 

Y-i-a^jm (is) 

Here, p is the area density of particles; and admittance Y is normalized to the admittance 
of free space. The 5(0) and S'(180) are the forward and backward scattering amplitude of 
a single particle, correspondingly, which are related to their electric and magnetic an 
susceptibilities as follows [42j H3] : 



5(0) = ik 3 (a E + a H ), 5(180) = ik 3 (a E - a H ) (14) 

The susceptibilities depend on particle concentration due to Lorentz field. Indeed, the 
susceptibility of a polarizable dipole in the array of identical dipoles is 

a = l^AH (15) 

Here cto is the susceptibility of an isolated particle (for a small isolated ideally conducting 
sphere, a E = r 3 ,an = — r 3 /2); a = (7rp) _1//2 is the Wigner-Seitz radius, and A is a local 
field factor [44] . To find refraction index and admittance of a planar ordered array of ideally 
conducting spheres, we substitute Eqjl5] into Eqs. T2|l3|14| and find 



x 2 (l + 2Ax 3 ) 

" 1+ (l-^3)(2 + Ax3) (16) 

Y ~ 1 + — (17) 

(1 - Ax 5 ){2 + Ax 3 ) 1 ' 



where x = r/a. Note that n, Y depend on density (through x and a). 
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To find refraction index and admittance of the planar disordered array of ideally con- 
ducting spheres we represent it by the ordered array with the spatially-varying density 
P — Po + $P- The Taylor expansion of Eqs T6fT7 with respect to p yields 

dn d 2 n So 2 

n( P )^n[ Po ) + -5 P+ —^... (18) 

^i^w+^Vt- (19) 



To find n, Y in the disordered state we use the same Eqs. 18fl9 and substitute Sp 2 by area 
density fluctuations. Then we average over the area S which scatters coherently. Its size is 
S = zX where z is the distance to the observation point and A is the wavelength. (When 
the averaging is performed over whole array, < 5p >= 0, but < 5p 2 >^ 0). 

The above approach reduces the effect of disorder on n, Y to the area density fluctuations. 
When the array of particles is ordered, this corresponds to compressional deformations. 
However, disorder introduces shear deformations as well. These deformations change local 
crystalline symmetry and this affects refraction index and admittance through the local field 
factor A ^QJ H3]. In principle, the effect of shear deformations may be treated along the 
same lines by considering Taylor expansion of Eqsfl6jT7] with respect to A. 



Electrodynamic parameters of the multilayer with in-plane disorder 

We represent our sample as a multilayer consisting of alternating layers of conducting 
spheres and air spacings between them. Wave transmission through multilayers is most 
easily accounted for by the matrix method (HI H5J. Here, each layer is characterized by two 
matrices: (i) the phase matrix representing the phase shift upon transmission through this 
layer, and (ii) the reflectivity matrix which is determined by the ratio of admittances at 
the reflecting interface and is independent of the layer thickness. Transmission through the 
multilayer is the product of all these matrices. We relate the effect of disorder on transmission 
through our multilayer to the area density fluctuations in the layers of conducting spheres. 
We assume that the lateral size of important fluctuations corresponds to the area that 
radiates coherently, as seen from the next interface, i.e. S ~ Xd z where d z is the unit cell 
period in the direction of propagation and A is the wavelength. In what follows we consider 
the effect of area density fluctuations in different frequency ranges. 
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Antigap 



The antigap is the narrow frequency range in the passband where the layer of spheres may 
be represented as a half- wavelength plate, i.e., the phase shift on propagation through this 
layer is a muliple of tt. Transmission through such layer is close to unity due to destructive 
interference of the waves reflected from both interfaces and is almost independent on the 
layer admittance. Transmission at the antigap frequency is, therefore, determined mostly by 



refraction index fluctuations across the layers (Eq, 18). The latter result in the phase shift 
fluctuations, if) = ipo + 5if). The power transmission through the multilayer decreases by the 
Debye- Waller factor, T disorder /T order = e~ 5 ^ 2 [37J. What is important here is the relative 
deviation of the phase shift when going from one layer to another. Hence, the averaging 
in the exponent is for different layers, i.e. for different realizations of the disordered state. 
Since if) = kond, where ko is the free space wavevector, then Sip = kodSn = kod^Sp. Hence, 
the ratio of transmittances for the iV-layer stack is 

ln^) = - N (kd^X W (20) 

J- order J antigap V "'P J 

Note that the ratio of transmittances is oc k$, i.e. increases with frequency. Therefore the 
effect of disorder is more pronounced in the higher-order passbands. 



Midgap 

In the midgap the reflected waves from all interfaces sum up in phase, hence the re- 
flectivity is maximal and transmission is minimal. For simplicity, we restrict ourselves to 
the first stopband and to the quarter-wavelength stack. Transmission through such stack is 
T f=a 1/Y 2N [16] where N is the number of layers. The average admittance of each layer in 
the disordered state is Ydisorder = Yorder + ^r^P 2 where the averaging is performed over the 
whole layer. This results in 

7^ \ ]\J M y 

i -'-disorder \ iy u I _ „ /ni\ 

^vzrj^-Yw 1 ' (21) 

Stopband edges 

To account for the effect of disorder on the transmission at gap edges we draw analogy 
to the optical properties of disordered semiconductors. Upon increasing disorder, the gap 
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in semiconductors becomes smeared and there appear exponential transmission tails at gap 
edges (Urbach tails). The frequency dependence of optical transmission through semicon- 
ductors and at the gap edge is [48J: 

dh\T 



df 



oc 5a\ ( 22 ) 



where 5a L is the lattice constant variation which can arise from structural or temperature 



fluctuations. We rewrite Eq.22 for a single frequency, and note that Scll/cll = 5p/2p. Then 



Eq.22 yields the ratio of the transmissions in the ordered and disordered states at any 



frequency corresponding to the gap edge as 



^^sardeA (x6p2 ^ 

J order J e d ge 



Scattering 



In addition to the frequency-selective effects considered in previous subsections, there is 
also scattering that decreases transmission at all frequencies [231 09] . These extinction losses 
are also proportional to density fluctuations [32]: 



T, 



In ^ sorder j oc5p 2 (24) 

order J scattering 



We observe that Eqs 20|2~T|22[24 - all yield exponential dependence of the transmission 
on density fluctuations, 



\n-^>=-B{f)^^> (25) 



T(f,H) _ 6?(f,H) 

Tord P 1 

where T ord is the transmission in the completely ordered state and B(f) is the frequency- 
dependent prefactor. Equation [25] assumes monotonous dependence of transmission on dis- 
order. This is valid only for the relatively weak disorder, while for the strong disorder this 
dependence can be nonmonotonous in certain frequency ranges [19]. 



COMPARISON TO EXPERIMENT 



Magnetic-field dependent order-disorder transition 

We explore magnetic field dependence of the in-plane disorder. We quantify disorder 
through the intensity of diffuse rings in the Fourier transform images (Figjl]). As it is well- 
known from the X-ray structure analysis, Idiff is directly related to the lattice constant 
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fluctuations, which are closely related to the density fluctuations, i.e. i^i// oc 5p 2 /p 2 . At 
Figs 8j9 we plot Idiff in dependence of magnetic field. Here, we set Idiff = 1 in the com- 



pletely disordered state, i.e. at H = 0. In the weak field the data can be well approximated 
by the linear dependence, Idiff = 1 — /3H as predicted by Eq{TT]for the amorphous state. 
At higher field, the data are fairly well approximated by the 1/H 2 dependence (see also 
Fig|9]), as predicted by Eqj9]for the polycrystalline state |60j. The crossover between these 
two regimes occurs at H = 50-60 Oe. This corresponds to the field when the Bragg peaks in 
the Fourier transform image are replaced by diffuse rings (Fig.l) and is suggestive of some 
kind of melting or glass transition [501 EH E2], |53j [M] • 

Mm-wave transmission through the regular crystal 

To verify relation between transmission and area density fluctuations, as predicted by 
Eq,25, we plot them together versus magnetic field (Figs, 8f9 ). Both Idiff (H) and T(H) 



dependences collapse on one curve. This means that T(H) is indeed proportional to /, 



and to the density fluctuations, as predicted by Eq. 25 This is the central result of our paper. 

The absolute value of the fluctuations in the completely disordered state may be estimated 
as follows. The Eq. 11 predicts that at H = 0, Sp 2 /p 2 ~ 1/Zs- Since Zs = pd z X, we 



substitute p = 0.17mm~ 2 , d z = 4.5mm, A = 6-15 mm, and find Zs ~ 5 — 11 depending on 
frequency. This corresponds to 5p 2 /p 2 = 0.1 — 0.2. 

Since our model accounts fairly well for the magnetic field dependence of the mm-wave 



transmission, in what follows we only estimate the prefactor B(f) in Eq.25 To this end we 



compare the ratio of transmittances in the completely ordered and in the completely disor- 



dered states and for ceratin frequencies. In particular, for antigap, we substitute into Eq.20 
N = 6, / =46 GHz, Jp^/p 2 = 0.2, a = 1.37 mm, x = 0.73, A = 0.5, n = 1.4, pdn/dp = 0.7 
and find lnTdi sorc fe r /T or< fe r =-24 dB. This should be compared to the experimentally observed 



value of -10-15 dB (Figj4j). For the midgap, we substitute into Eq. 21, iV = 6, / = 27 GHz, 
Jffijp 2 = 0.12, x = 0.73, A = 0.5, Y = 1.9, p 2 d 2 Y/dp 2 = 0.4 and find lnT disorder ./T or(ier .=-2.5 
dB. This should be compared to the observed value of ~-5 dB (Figs 3p ). 



We conclude that our model qualitatively predicts field dependence of the transmission 
through disordered photonic crystal. To achieve better quantitative agreement, there is need 
in a more advanced model of the kind suggested by Ref.[55j. This model should go further 
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than Eqs, T2fT3~|l6||17 which assume low concentration of scatterers and Rayleigh scattering 
regime. The model should also account for the quadrupole interactions between magnetized 
particles. 



Mm-wave transmission through the crystal with a planar defect 

The microwave transmission through the crystal with a planar defect can be treated 
following the same lines. We consider the photonic crystal split on two halves, as two 
photonic crystals coupled through the Fabry-Perot resonator. Transmission through this 
device is [31]: 

where A is the spacing between the mirrors and S12 and Su are complex transmission and 
the reflection coefficients for each of the halves (note that Sn and S12 are interrelated and 
both are affected by disorder). The resonant transmission through such a crystal is possible 



rr JiZ a ~, (26) 



provided the denominator of Eq.26 is small, i.e. 1 — |Sii| 2 <C 1. Even small disorder strongly 
affects this condition, hence the transmission through the defect drops dramatically in the 
disordered state (Fig.??). 



DISCUSSION AND CONCLUSIONS 



Continuous variation of electromagnetic properties of our samples when external mag- 
netic field goes to zero closely resembles the change of optical properties of liquids upon 
approaching the critical point. Indeed, the optical transmission at the critical point is very 
small due to enhanced density fluctuations. According to Einstein's theory, the density fluc- 
tuations inversely depend on compressibility which goes to zero at the critical point [IT] . 
The compressibility of our array of magnetizable spheres, K^ 1 ~ 1/H 2 (EqjSJ), also diverges 
when H — > 0. The density fluctuations, as it is shown by Eq|9j increase accordingly. 

Our experiments bear some resemblance to the studies of electromagnetic wave propa- 
gation through ferrofluids [56j EZ] • Indeed, in the absence of magnetic field the particles in 
ferrofluid are in the disordered state, while in the presence of magnetic field they self-assemble 
into chains. This results in field-induced anisotropy which strongly affects the polarization 
of the electromagnetic wave propagating through such a media. However, magnetic field 

16 



does not induce the long-range order in ferrofluids, hence the field-induced anisotropy is 
frequency-independent. This is very different from our present experiment where magnetic 
field drives the particle array into crystalline state with a long-range order, and this has a 
strong frequency-selective effect on the microwave transmission. However, two-dimensional 
ferrofluid layers j58j [59] bear strong resemblance to our system. 

In summary, we demonstrate a metallo-dielectric photonic crystal exhibiting stopband 
in the mm-wave range. Mm-wave transmission through this crystal can be controlled by 
external magnetic field through magnetic field- induced order-disorder transition. We develop 
a physical model which describes our experimental data fairly well. Our concepts may be 
useful in interpretation of the effect of disorder in photonic bandgap materials. Our results 
can be useful for the fabrication of tunable planar photonic crystals, based on surface waves, 
in particular surface plasmons. If the surface is covered with the liquid layer containing 
movable magnetic particles, (especially when this layer represents a photonic crystal with a 
tunable lattice constant) the propagation of the surface waves can be effectively monitored. 
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raeli Ministry of Science and Technology. We are grateful to I. Felner for magnetization 
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FIG. 1: Particle arrangement in a container at different values of magnetic field. The container is 
made of a 0.67 mm thin plexiglas plate with the hexagonally-shaped walls. The side of a hexagon 
is 30 mm and there are 397 steel spheres of 2 mm diameter. The left panel shows grey scale 
images of the particle configuration as obtained by the CCD camera. [The illumination was from 
above, therefore the actual diameter of the spheres exceeds the diameter of the small white circles 
in the real space images. In reality, the spheres in Fig. lc touch one another.] Note gradual 
transition from the ordered to disordered state upon decreasing magnetic field. The right panel 
shows corresponding Fourier transform image. Note gradual disappearance of sharp peaks and 
appearance of diffuse rings upon increasing disorder. 
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FIG. 2: Measurement Setup. A stack of 6-10 layers with steel spheres is mounted inside Helmholtz 
coils which produce magnetic field perpendicular to the layers. The mm- wave transmission through 
the stack is measured using standard gain microwave horns connected to HP 8510C Vector Network 
Analyzer. The antennas are equipped with collimating teflon lenses. 
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FIG. 3: Mm-wave transmission through the six-layer stack for different values of magnetic field. 
The interlayer separation is d = 3.5 mm. Note stopband at 25-44 GHz and at H =290 Oe (ordered 
state) and its disappearance in the disordered state (H =0). 
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FIG. 4: Mm-wave transmission for the similar stack but with increased interlayer distance (d = 
4.5 mm). In the ordered state {H =290 Oe) note first stopband at 21-36 GHz and the edge of 
the second stopband at 46 GHz. Note smearing of the stopband in the disordered state and broad 
transmission peak at 39 GHz (antigap). 
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FIG. 5: Magnetic field dependence of the mm-wave transmission through the sample of Fig. 4 and 
at fixed frequencies. / =34.7GHz corresponds to the high-frequency edge of the first stopband, 
/ =39.1 GHz corresponds to the antigap, and / =46.3 GHz corresponds to the low-frequency edge 
of the second stopband. 
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FIG. 6: Mm-wave transmission for the same stack as in Figj3] but with planar defect (the stack 
shown in FigJTjwas split in two parts separated by 6 mm). Note a sharp peak at 34 GHz inside 
the stopband which corresponds to the resonant transmission. 
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FIG. 7: Dispersion relations for the ten-layer stack (331 sphere in each layer, interlayer distance 
is 4.3 mm). Filled symbols show Re(/c) and open symbols show Im(k), where k is the wavevector. 
The dashed line shows ko- the wavevector in the uniform medium with the same refraction index 
as our photonic crystal, (a) The ordered state, H =160 Oe. (b) The disordered state, H =0. Note 
the gap between 24 GHz and 36 GHz in the ordered state and its absence in the disordered state. 
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FIG. 8: Mm-wave transmission at 46.3 GHz for the sample of Fig. 5 (filled circles) and the intensity 
of diffuse background (open rhombs), Idiffi estimated from the Fourier-transform images of Fig.l. 
Straight solid line shows linear dependence dependence as predicted by Eq{TT]for the glassy state; 
the thin curved line shows 1/H 2 dependence as predicted by Eqj9]for the polycrystalline state. The 
vertical line delineates the strong field region (where Fourier transform imaqe of Fig. 1 shows sharp 
Bragg peaks, indicating on crystalline state), from the weak field region (here Fourier transform 
imaqe of Figjl] shows sharp diffuse rings indicajsgig on amorphous or glassy state.) 




FIG. 9: A different representation of the high-field data of Fig|8] Filled circles show mm-wave 
transmission at 46.3 GHz for the sample of Fig. 5, open rhombs show the magnitude of area density 
fluctuations Idiff- Straight solid line shows linear dependence on 1/H 2 as predicted by Eq|9j 
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